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accelerated universe. The role of physical parameters have also been shown graphically during 
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I. INTRODUCTION 

Cosmic acceleration is on of the most challenging observation of the cosmology 0, 0. The reason for 
this accelerating universe is termed as dark energy which dominates the universe (70% of the universe) 
having large negative pressure and violates the strong energy condition The most popular candidate 

of dark energy is the cosmological constant A Q whose EoS parameter is w = — 1. This kind of dark 
energy shows that the universe will accelerate forever. Another kind of dark energy dubbed as quintessence 
(w < —1/3) explains that the acceleration will replaced by deceleration in far future. And for phantom 
energy (w < —1), acceleration will change to super acceleration which will eventually destroy every stable 
gravitational structure ■ There are also other candidates of dark energy like Chaplygin gas [1] , modified 
Chaplygin gas (MCG) Tachyonic field fiof . DBI-essence 11 1, K-essence [l2| and so on. 



There are numerous works done on dark energy on the theory of Einstein's classical general relativity 
(GR). But, most physicists deemed that the gravity should be quantized. Loop quantum gravity (LQG) 
is an outstanding effort to describe the quantum effect of our universe. In this theory classical space time 
continuum is replaced by discrete quantum geometry. Now a days several cosmological (interacting dark 
energy model) models are studied in the frame work of LQC. Wu and Zhang [l3[ studied the cosmological 



evolution in LQC for the quintessence model. Chen et al [14j provided the parameter space for the existence 
of the accelerated scaling attractor in LQC with more general interacting term. When the Modified 
Chaplying Gas coupled to dark matter in the universe is described in the frame work LQC by Debnath et 
al |15| who resolved the famous cosmic coincidence problem in modern cosmology. 

There is another modification on gravity (Brane-gravity) which also exhibits the acceleration of the 
present day universe. It is proposed that our universe is a 3-brane embedded in a four dimensional space. 
An important ingredient of the brane world scenario is that the standard matter particles and forces are 
confined on the 3-brane and the only communication between the brane and bulk is through gravitational 
interaction (Le., gravity can freely propagate in all dimensions) or some other dilatonic matter. In the 
review [!&Jl9| there is different applications with special attention to cosmology in Brane-gravity. In this 
work we consider the two most popular brane models, namely DGP and RS II branes. 

Regarding cosmological acceleration, dark energy with energy density of scalar field act subdominant 
during radiation and dark matter eras and acts dominant at late times. Dynamical system theory has been 
applied with great success in cosmology and astrophysics within the context of general relativity. This 
theory are used to describe the behaviour of complex dynamical systems usually by constructing differential 
equations. This theory deals with a long term qualitative behaviour of the formed first order differential 
equations. It does not concentrate to find the precise solutions of the system but provide answers like 
whether the system is stable for long time and whether the stability depends on the initial conditions. 
Besides the other scientific fields this theory is now become widely useful in the research of cosmol ogy . In 
the construction of different dark energy model cosmological scaling solutions work significant role |20h23| . 
Tsujikawa et al 0, [24| proved that scaling solution exists for coupled dark energy whenever they restrict 
to the form of the field Lagrangian p(X, </>) = Xg(Xe x< ^) where X = —g^ v d^<j)d u (j)/2 and g is any function 
of Xe x ^. In reference (2||, they also considered the interacting model in these Lagrangian form and studied 
the stability of fixed points for several different dark energy models for ordinary (phantom) field, dilatonic 
ghost condensate and (phantom) tachyon. Our main aim of this work is to examine the nature of the 
different physical parameters for the universe around the stable critical points in LQC and two brane world 
models (DGP and RS II) in presence of DBI-esssence type dark energy along with dark matter with suitable 
interaction term. With the evolution of the universe we find the effective state parameter w e ff, Critical 
densities for dark energy (f2^) and for dark matter (f2 m ) and examine future dominance nature of kinetic 
energy and potential energy. 

In this work, we have considered the field Lagrangian p(X, 4>) — Xg(Xe x ^) and studied the dark energy 
model in different gravity theories like (i) LQC, (ii) DGP Brane-world, (iii) RS II Branc-world. We also 
derive the critical point of the dynamical system in different gravity and analyze the stability. Also we 
do the numerical simulation for LQC, DGP-Brane-world and RS II Brane-world models. Some fruitful 
conclusions are drawn in section V. 
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II. BASIC EQUATIONS IN DBI-ESSENCE 

The action of the Dirac-Born-Infeld (DBI) scalar field (f> can be written as (choosing 8irG = c = 1) [2f 



(i) 



where V(<fi) is the self-interacting potential and T(cj>) is the warped brane tension. The kinetic term of 
the above action is non-canonical. Physically, this originates from the fact that the action of the system 
is proportional to the volume traced out by the brane during its motion. This volume is given by the 
square-root of the induced metric which automatically leads to a DBI kinetic term. 

From the above action, it is easy to determine the energy density and pressure of the DBI-essence scalar 
field which are respectively given by 

^ = ( 7 -l)T(0) + y(0 (2) 

and 

P4 , = ^^T(4>) - V{<t>) (3) 
7 

where 7 is given by 

^ (4) 
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From above expression, we observe that T(<f>) > <fi 2 and thus 7 > 1. We consider a spatially flat Friedmann- 
Lemaitre- Robertson- Walker (FLRW) Universe containing a perfect fluid and a scalar field (j>. Assuming that 
there is an interaction between scalar field (dark energy) and the perfect fluid (dark matter) , so they are not 
separately conserved. The energy balance equations for the interacting dark energy and dark matter can be 
expressed as [HI 

P0 + 3H(1 + w^)p0 = ~Qp m 4> (5) 

and 

f> m + 3i?(l + W m )p m = Qpm4> (6) 

where p m is the energy density of the dark matter, w m is the EoS parameter for the dark matter, H = — 
Hubble parameter, a is the scale factor and Q > is the coupling between dark energy (DBI-essence) and 
the dark matter. 

We define the fractional density of dark energy and dark matter, = t-^j- and f2 m = ■Sffr ■ 

To get stable attractor solution we must have 7 = constant, — Xg(YY is the scalar field pressure 
density where Y = Xe X( ^, X = —g^d^dv^jl with g is any function of Y [23|, [Hj]. For DBI-essence we 
choose 

T(0) = V{4>) = V e^ (7) 

Then the pressure p^, and the energy density p^ can be written in the form 

= Xg(Y) 

P<p = 2X^ - P4> = X[g{Y) + 2Yg'(Y)} (8) 
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whenever we choose 



where ' denotes the derivative with respect to Y. The total cosmic energy density p = p^, + p m satisfies the 
conservation equation p + 3H(p + p) = 0, where p = p<p + p m - 



III. EVALUATION OF DYNAMICAL SYSTEM IN LQC 

The modified Friedmann equation for LQC is given by [l3|, Q, H3| ■ 



H 2 = P - 



(10) 



Here p c = \ // 3n 2 rj' i G 2 h is the critical loop quantum density and rj is the dimensionless Barbero-Immirzi 
parameter. It should be noted that for our LQC model, p < p c . 



Consequently we obtain the modified Raychaudhuri equation (using the conservation law) 

H = -\{p + p){l-2^) 
2 p c 



(11) 



We introduce the following dimensionless quantities 



V6H' 



y = 



,-X<t>/2 



(12) 



We see that y and z(< 1) must be non-negative, but x may or may not be positive depends on the nature 
of <j>. Substituting the expressions of x, y and z in equations (5), (6) and (11), we obtain the first order 
differential equations in the form of autonomous system as follows: 



dx 3x 



A(l - w m )x 



(l+Lu m ){l + V y 2 (z-l)} 



z- 1 



(1 - 2z) + 



V6Qx 2 {1 + Ax 2 (z - 1) + V y 2 {z - 1)} V6XV y 2 (3x 2 - 4y 3 )(l - z) 



dy_ _3y 
dN ~ 2 



2Ax 2 {z-\) 2Ax 2 (z-\) 
(l+w m ){l + V y 2 (z-l)} 



A(l - w m )x - 



dz 
dN 

1 dH 

HdN 



= -3 
3 

~ _ 2 



A(l - w m )x 2 - 
A(l - w m )x 2 - 



z- 1 

(l + w m ){l + V y 2 (z- 1)} 
z- 1 

(l + w m ){l + V y 2 (z-l)} 
z-1 



(1 - 2z) ^Xxy 



(l-z)z 
(1 - 2z) 



(13) 
(14) 

(15) 

(16) 
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where A = ry 2_ 1 > N — In a being the number of e- folds and a being the scale factor. In terms of the new 
variables x, y, z we obtain the following physical parameters 



V y 2 , 



W eff 



Ax 2 - V y 2 
Ax 2 + V Q y 2 

w m + (Ax 2 - V y 2 )(l - z) 

P4, + p m l + Ax 2 (l-z) + V y 2 (l-z) 



^0 = Ax 2 H 

* + Pn 



(17) 
(18) 
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x(N) y(N) *(N)| O 1 J 

FIG.l FIG.2 



FIG. 1: Evaluation of x, y, z with respect to N in LQC model taking 7 = 1.2, Q — 0.05, Vo = 8, A = 0.5 and 

W m = 0.01. 

FIG. 2: Evaluation of O^,, f2 m and w e // with respect to iV in LQC model taking 7 = 0.8, Q = 0.05, Vo = 7.6, A = 0.5 
and w m = 0.01. 

The fraction density function as = jj/p and fi TO = ^Sj- satisfying 

+ n m + Qlqc — 1 

where, Qlqc = p _f p is the density parameter due to the effect of LQC. Since p < p c , so ^lqc < in our 
case. 

The new variables [x, y, z) have been drawn in figure 1 with respect to N = In a and seen that all are 
positive oriented due to the expansion of the universe. Also fi^, f2 m and w e ff have been drawn in figure 2. 
f2 m shows the lower value (< 1) and f2^ shows the value > 1 in evolution, so fi^ gets higher value than VL m . 
So in late stage, the dark energy (DBI-essence) dominates over dark matter. Also w e ff gives the negative 
value less than —0.5 which shows the dark energy dominated phase of the universe. 



1. Critical points: 



The critical points can be obtained by setting 
following table. 



dx 
JN 



0, & = and ft 



and are presented in the 
Table 1: The critical points (x c , y c , z c ) and the corresponding values of the density parameter VL^. 



No. x c 



Vc 



fid 



1 



\ m ) A(w m -1) 



(iv) 



2Q 
















2Q 2 



3A(1— wU 



3A(l-w m ) 2 
3A{l+w m ) 



2Q 2 



From the Table 1, we see that the components of y c is equal to zero for above four critical points. The value 
of fi^ = 1 for the critical points given in (i) and (ii). These provide the accelerated phase of the universe. 
Similar nature happen for other two critical points (iii) and (iv), but these depend on the interaction term 
Q, A and w m . 
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2. Stability of the model: 



Now the stability around the critical points can by determined by the sign of the corresponding 
eigen values. If the eigen values corresponding to the critical point are all negative, the critical points are 
stable node, otherwise unstable. The eigen values for the above critical points are obtained as in the following: 

Table 2: The eigen values corresponding to the critical points (x c , y c , z c ). 



No: Valuel Value2 Value3 



(i) 

(ii) 
(iii) 
(iv) 



V6Q 



3«v, 



3 _)_ 3^77 

' 2 ' 2 



Q 2 

A(l-w m ) 



-^2! 3(1 + Wm ) 



■l{l + R) 



A(l-w m ) 



K-i + h) 



n _ / 3A 
V 2A 



3 + 



3A(1— <J+2Q(Q+A) 
2A(-l+w m ) 

3(l+w m )\ 
2Q 



wlu , rei R = J 6 A(-l +Wm )(l + ^-QH3 +iWm ) 



From the Table 2, we observe in the following: 

(a) One eigen value for the critical point (i) is positive, since 3(1 — w m ) + -y=2 > 0, so around this critical 
point system is not stable. 

(b) One eigen value for the critical point (ii) is positive, since, 3 + \J~^ > 0, so around this critical point 
system is not stable. 

(c) If 7 2 < 1 and Q(Q + A) > 3A(1 ~ W "' ) , the eigen values of the critical point (iii) are all negative and 
hence the system is stable (node). 

(d) All eigen values of the critical point (iv) are negative if A < and R < 1, so the system may be stable 
otherwise the system will be unstable. 



IV. BRANE WORLD 



A. Basic equations in DGP Brane model 



An effectual model of brane-gravity is the Dvali-Gabadadze-Porrati (DGP) braneworld model [28, 29] that 
represents our 4-dimensional universe to a FRW brane embedded in a 5-dimensional Minkowski bulk. It 
explains the origin of DE as the gravity on the brane escaping to the bulk at large scale. On the 4-dimensional 
brane the action of gravity is proportional to Ml . That action is proportional to the corresponding quantity 
in 5-dimensions in the bulk. The modified Friedmann equation in DGP brane model considering flat, 
homogeneous and isotropic brane is given by 



(19) 
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M 

where p is the cosmic fluid energy density, H = |, Hubble parameter and r c — ^yz- is the crossover scale 
which resolve the transition from 4D to 5D behavior and e = ±1. Corresponding to e = +1 the we have 
standard DGP(+) model which is self accelerating model without any form of DE, and effective w is always 
non-phantom. However for e = —1, we have DGP(— ) model which does not self accelerate but requires DE 
on the brane. Using (18), the modified Raychaudhuri equation becomes (choosing 8ttG = c = 1) 



2H-e- )H = -H(p + p) 

T ' c. 



(20) 



1. Dynamical system 

To get dynamical analysis of our DGP brane world model of the universe, we define the following dimen- 
sionless quantity 



V6H' 



y = 



y/ZH 



z = 2H 



(21) 



We see that y must be non-negative, but x and z may or may not be positive. When <f> is increasing, x 
must be positive and <j> decreases implies x is negative. Also e = — 1 implies z > 0, but for e = +1, the 
value of z may or may not be positive. Now we introduce the fraction density parameters as fl^ = and 
= j^T satisfying 



(22) 



where, Qdgp = 7-77 ^ s ^ e density parameter due to the effect of DGP brane world with the physical 
parameters 



0,4, = Ax 2 + V y 2 , wj, = 



Ax 2 - V y 2 



Weff = 



Ax 2 + V y 2 

+ p m _ 4r 2 c (Ax 2 - V y 2 ) + w m [-1 - 4r c 2 (-l + Ax 2 + V y 2 ) - 4r c e + e 2 )} 



P4> + Pm (2r c + e) 2 - 1 

Using all this and defining N = In a (the number of e-folds) we get the system of equations as follows: 



(23) 
(24) 



dx „ 3x ( e \ / t 9 
— = -3a; + — (z + — 2Ax 2 
dN 2z V r c J V 



(1 + w m ){l + 4r c 2 (-l + Ax 2 + V Q y 2 ) + 4r c e - e 2 } 
4r 2 



V6[Q{1 + 4r e 2 (-l + Ax 2 + V y 2 ) + 4r c e - e 2 } + 4U Ar c 2 y 2 ] 

+ 8M (25) 



dy 
dN 



I {z + ±) (2Ax 2 - (^^^r 2 (-l + Ax 2 + V y 2 ) + 4r c e-e 2 ^ + ^ 



(26) 



dz 3 / e 

diV = ~2~z \ Z + 7 r 



[2Ax 2 



(1 + w m ){\ + 4r 2 (-l + Ax 2 + U y 2 ) + 4r c e - e 2 } 
4r 2 



(27) 



1 dH 
HdN 



3 / 



2Ax 



2 (1 + w m ){l + 4r e 2 (-l + Ax 2 + V y 2 ) + 4r c e - e 2 } 

4r 2 



(28) 



The new variables (x,y,z) has been drawn in figures 3 and 5 with respect to N = In a for DGP(+) 
and DGP(— ) models respectively. In all the cases, x has shown to be negative i.e., the DBI scalar field 4> 
decreases during expansion, but y and z keeps positive sign. The effective EoS parameter w e ff and the 
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density parameters fi^, f2 m are shown in figures 4 and 6 for DGP(+) and DGP(— ) models respectively. 
During expansion, the £7^ increases and fi m decreases which show the dark energy dominates at late times. 
Also w e ff decreases from some —0.4 to —1 for both models, which also shows the dark energy dominated 
phase of the universe. 




y(If) z(N) 



-0.6- 
-0.8 









£2 
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V 




/ N 




/ \ 




/ \ 




/ - v 






e - 1 




\ 





FIG.3 



FIG.4 



FIG. 3: Evaluation of x, y, z with respect to N in DGP (+) Brane model (e = l)taking 7 = 0.6, Q = 0.05, Vq = 0..' 
A = 0.5, w m = 0.01 and r c = 1000. 



FIG. 4: Evaluation of f2^, Q m and w B ff with respect to N in DGP (+) Brane model (e = 1) taking 7 = 0.6, Q = 0.05, 
V = 0.8, A = 0.5, w m = 0.01 and r c = 1000. 




FIG. 5: Evaluation of x, y, z with respect to TV in DGP (-) Brane model (e = — l)taking 7 = 0.6, Q — 0.05, Vb = 0.8, 
A = 0.5, w m = 0.01 and r c = 1000. 

FIG. 6: Evaluation of Q. m and w e ff with respect to N in DGP (-) Brane model (e = — l)taking 7 = 0.6, Q = 0.05, 
Vb = 0.8, A = 0.5, w m = 0.01 and r c = 1000. 
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2. Critical Points: 



The critical points can be obtained by setting = 0, 4& — and jjjf — 0. The possible critical points 
(x c , y ci z c ) and the corresponding values of of our DGP model are given by 



(i\fn / 0{-l+(e-2r^T e\ o _ Q{-l+(e-2r c ) 2 } 
V) y 4V (Q+\)r 2 ' r c J' >,l <f>- 4(Q+\)r 2 

r ..s / ^Ar 2 + y /Ar 2 l6Ar 2 +Q 2 {-l + (e-2r7) 2 }] n e \ _ ( vgArg+ygg [&Ar 2 +Q 2 {- l + (e-2r c ^}J) ° 

Wl 2AQr 2 ' U > — r c ) ' ^4> — 4AQ 2 r% 

V6Ar 2 - v /Ar 2 [6Ar 2 +Q 2 {-l + (c-2r c ') 2 T] £ \ _ (- ygA)j+ygg [6Arg+Q^ {- l + (e~2r c nj) ' 



\ 2AQr 2 r c )' QD< P 4AQ 2 r* 

The value of fi^ > 1 or < 1 for the critical points given in (i) to (iii), depends on the values of the r c , A 
and the interaction term Q. 

3. Stability of the model: 

Now the stability around the critical points can by determined by the sign of the corresponding 
eigen values. If the eigen values corresponding to the critical point are all negative, the critical points are 
stable node, otherwise unstable. The eigen values for the above critical points are obtained as in the following: 

Table 3: The eigen values corresponding to the critical points (x c , y c , z c ). 



NO: 


Valucl 


Value2 


Value3 















3Ar 2 + ^/9A 2 ri-3AQr 2 {-l + (-2r c +e) 2 }\ 


-3Ar 2 + v /9A 2 r4-3AQr 2 .{-l + (-2r c +e) 2 }\ 


(i) 


2Ar 2 


2Ar 2 


(ii) 





v /3[Ar 2 {6Ar 2 +Q 2 {-l+(-2r c +e) 2 }}] 


{QAr 2 + ^/Q[Ar 2 {QArl+Q 2 {-l+(-2r c +e) 2 }}])\ 


V2Ar 2 


4Ar 2 


(iii) 





v /3[Ar 2 {6Ar 2 + Q 2 {-l+(-2r c +e) 2 }}] 


(-6Ar 2 + y?6lAr 2 {6Ar 2 +Q 2 {-l + (-2r c +e) 2 }}])\ 


V2Ar 2 


4Ar 2 



From Table 3, we see that one eigen value for all three critical points is zero. Hence the dynamical system 
is unstable around all critical points. 

B. Basic Equations in RS II Brane World 

Randall and Sundrum [3"ill32T[ elucidate the higher dimensional scenario by introducing a bulk-brane model 
dubbed as RS II brane model. They proposed that we live in a four dimensional world (called 3-brane, a 
domain wall) which is embedded in a 5D space time (bulk). All matter fields are confined in the brane and 
gravity can only propagate in the bulk. In RS II Brane world the modified Einstein equations in flat universe 
are 



■AH 2 = A 4 + n\p + ^p 2 + " W (29) 

2H + 3ff 2 = A 4 - n\p - ^pp - ^L 2 - * W (3Q) 

zAi l\\ 
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Here K4 and A4 are respectively 4D gravitational constant and effective 4D cosmological constant. The 
dark radiation U satisfies the relation 



U + 4HU = 



(31) 



1. Dynamical system 



Here we introduce the new variables 



<p _ e "v" p 



(32) 



Also we introduce the fraction density function as toj, = ^7 and to m = -^ffj satisfying 

n,p + to m + ansa = 1 (33) 

where, tonsil = 1 sttx \ is the density parameter due to the effect of RS II branc world with, 

Ax 2 - V y 2 



to^ = Ax + V y , = 



Ax 2 + V Q y 2 



W eff = Kl{l + z) 



A(l - w m )x 2 - V (l + w m )y 2 + 



K 2 4 {l + Z) 



(34) 

(35) 



In absence of the cosmological constant and dark radiation (A4 =14 = 0) the above equations reduce to 
the dynamical system of equations as follows 



dx q ^ 2 
— = -3x+- K4 x 



(Ax 2 + V a y 2 ) z + 2 Ax 2 (I + z) + {1 + w m + (2 + w m )z} (-Ax 2 - V y 2 + 
V3 [-Q (Ax 2 + Voy 2 ) + - V \y 



1 



nl(l + z) 



V2A 



(36) 



dy = y 
dN _ 2 



l + z 



+ SAk\(\ - w m )x 2 (l + z) - 3K 2 V Q y 2 (l + z) - 3iu m {-l + n\V a y 2 (\ + z)} - v^Aa 



(37) 



dz 
dN 



-3kIz(1 + z) 



Ax 2 (I - w m ) + (1 + w„ 



K 2 (Z + 1) 



v y 



(38) 



1 dH 
HdN 



7^4 



(Ax 2 + V y 2 ) z + 2Ax 2 (l + z) + {1 + w m + (2 + w m )z} (-Ax 2 - V Q y 2 



1 



k 2 4 (1 + z) 



(39) 



The new variables (x, y, z) has been drawn in figure 7 with respect to N = In a for RS II model. We see 
that x, y, z are shown to be positive throughout the evolution. The effective EoS parameter w e ff and the 
density parameters to^, to m are shown in figure 8. During expansion, the to^ increases and to m decreases 
which show the dark energy dominates at late times. Also w e ff decreases from some value 0.1 to negative 
value < — 1 which also show the dark energy dominated with phantom phase of the universe. 
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FIG. 7: Evaluation of x, y, z with respect to N in RS II Brane model taking 7 = 0.6, Q = 0.05, Vo = 6.5, A = 1.5, 
Wm = 0.01 and ka = 1. 

FIG. 8: Evaluation of fi^,, Sl m and w e // with respect to N in RS II Brane model taking 7 = 0.4, Q = 0.2, Vo = 6.5, 
A = 1.5, w m = 0.01 and K4 = 1. 



2. Critical points: 

The critical points can be obtained by setting = 0, -fft = and ^ = 0. The feasible critical points 
(x c ,y c , z c ) and the corresponding values of il^ of our RS II model are given in the following: 

Table 4: The critical points (x c , y c , z c ) and the corresponding values of the density parameter VL^. 



No. x c 



z c 



(1) 



,..s V2Q+y / 2Q 2 -3AK 2 (w 2 n -lj 
V3A K 2 (w m -l) 



V6Q-j6Q 2 -9A K , 2 (w 2 n -l) 
V 11 ' 3AKi( Wm -l) U 



6Anl(~l+w 2 Tl )+2Q(2Q+^/4Q 2 -6AK,l(w 2 n -l)) 
3Ak 2 (w 2 -1) 



-6Ak,1(-1+w? 11 )+2Q(2Q+^/4Q 2 -6Ak 2 (w^-1)) 
SAk'Hw^-I) 



2Q 2 



3AK,\(w m -l) 2 



(V2Q+yf2Q 2 -3AK 2 (w 2 m -l)) ' 
3AK,\{w m -\) 2 



[--y2Q+^/2Q 2 -3AK 2 (wl l -l)j 
3Anl(w m -l) 2 



From the Table 4, we see that the components of y c are equal to zero for above three critical points. The 
value of f2^ > 1 or < 1 for the critical points given in (i) to (iii), depends on the values of the w m , A and 
the interaction term Q. 



3. Stability of the model: 



Now the stability around the critical points can by determined by the sign of the corresponding eigen 
values. If the eigen values corresponding to the critical point are all negative, the critical points are stable 
node, otherwise unstable. The eigen values for the first critical point (i) are obtained as in the following. 
The eigen values of critical points (ii) and (iii) are very difficult to obtain, so we have not considered that 
critical points here. 
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Table 5: The eigen values corresponding to the first critical point (x c , y c , z c ). 



NO: Valuel Value2 Value3 



m 2Q 2 o ( , , \ Q 2 3fi ,„ \ 3A^(-l+^, 2 „)-2Q(Q+A) 

W AKi(w m -l) °l i + ^mj Are^(tu m -1) 2 V 1 "'mj 2A K 2 (-l+u> m ) 



At the critical point (i), the three eigen values cannot be negative simultaneously, since Q and w m are 
small quantities. So the dynamical system is unstable in this case. At the critical points (ii) and (iii) we can 
not study the stability analysis. 



V. DISCUSSIONS 



In this work, we have studied homogeneous isotropic FRW model having dynamical dark energy DBI- 
essence with scalar field in presence of perfect fluid having barotropic equation of state (i.e., p m = w m p m ). 
The existence of cosmological scaling solutions restricts the Lagrangian of the scalar field <f>. The stable 
attractor solution can be found only for 7 = constant. We have chosen the potential function for DBI-essence 
as V(<j>) = V e Xc t>. Choosing p = Xg{Xe x *), where X = -g^ v 8^6^/2 with g is any function of Xe Xct > 
and defining some suitable transformations, we have constructed the dynamical system in different gravity 
theories like (i) Loop Quantum Cosmology (LQC), (ii) DGP Brane World and (iii) RS-II Brane World. 
For all gravity models, VL m gradually decreases to a small positive value and VL^ gradually increases to a 
value near about 1. That means, DBI dark energy dominates over dark matter in late times. Also from the 
figures of w e ff, we see that w e ff keeps negative sign in late times. For LQC model, w e ff lies between —0.5 
and —1, which is the dark energy dominated phase. For DGP model, w e ff lies between —0.4 and —1. So 
for LQC and DGP models of the universe, the DBI dark energy valid only for quintessence era, they can 
not generate phantom era. Also in RS II model, w e ff < —0.1 and decreases to —1 upto certain stage of 
time and after that stage w e ft becomes less than —1. So in RS II model, the DBI dark energy valid for 
quintessence era and phantom era in late times. We have found some critical points and investigated the 
stability of this dynamical system around the critical points for three gravity models and inve stig ated the 
scalar field dominated attractor solution in support of accelerated universe. Gumjudpai et al [25[ in their 
work considered the dynamical system analysis for phantom field, tachyonic field and dilaton models of 
dark energy. They analyzed different dark energy model of scalar field coupled with barotropic perfect fluid 
and depicts that scaling solution is stable if the state parameter > — 1 and the scalar field dominated 
solution becomes unstable. The fixed points are always classically stable for a phantom field, implying 
that the universe is eventually dominated by the energy density of a scalar field if phantom is responsible 
for dark energy. Therefore in this case the final attractor is either a scaling solution with constant 
satisfying < O^, < 1 or a scalar-field dominant solution with 0,^ = 1, For our LQC model, four critical 
points have been found, in which only two critical points may be stable node while all other two critical 
points are unstable. For DGP model, three critical points have been found but they are all unstable. Also 
for RS II model, the calculated critical point is also unstable. An attractor scaling is established by Martin 
and Yamaguchi [ll[ after considering a dark energy model with DBI field. Recent times, the model of 
interacting dark energy has been explored in the framework of loop quantum cosmology (LQC). On that 
framework an interacting MCG with dark matter has been studied by constructing a dynamical system and 
depicts a scaling attractor solution which resolve the cosmic coincidence problem in modern Cosmology [l5| . 
A dynamical system is explored in DGP, RSII Brane world separately with suitable interacting dark energy 
coupled with dark matter model [33[ and investigated that the universe in both scenarios follow the power 
law form of expansion around the critical point. So in conclusion, DBI-essence plays an important role of 
dark energy for FRW model of the universe in loop quantum cosmology, which drives the acceleration of 
the universe. 
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